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Comment on “One-Dimensional Disordered 
Bosonic Hubbard Model: A Density-Matrix 
Renormalization Group Study” 

In a recent letter Pai et al. [ilj have applied the density- 
matrix renormalization group (DMRG) rnethod to the 
study of quantum phase transitions in interacting disor- 
dered Bose system. They claim to obtain an accurate 
phase diagram for the Hamiltonian 
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containing hopping term, on-site repulsion and random 
potential /q uniformly distributed between — A and A. 
The authors proved that the rnethod works remarkably 
well in pure system and allows precise calculation of the 
phase transition point between the superfluid (SF) and 
Mott-insulator (MI) phases. This was considered as a suf- 
ficient ground to apply DMRG rnethod to the disordered 
system. In this comment we agrue that DMRG com- 
pletely fails in predicting phase diagram in disordered 
case. At considerable disorder it has wrong shape and is 
completely incorrect quantitatively, as proved by quan- 
tum Monte Carlo simulations. For weak disorder the 
phase diagram obtained in Ref. [|lj contradicts theoreti- 
cal predictions (and an exact theorem!). 



The phase diagram for the Hamiltonian (1). The error 
bars are within the points size and originate mostly from 
the finite-size corrections (the system size used in this 
calculations was L = 100 as in Ref. [dj]). BG phase at 
the point maked by the cross was verified for the system 
with L = 1000 and t/T = 200. Dashed line indicates the 
results obtained by the DMRG rnethod [Q . 

Recent developments in the quantum Monte Carlo 
(QMC) algorithms allow one to calculate directly the su- 
perfluid stiffness, compressibility and the dependence of 
the number of particles on the Chemical potential in sys- 
tems as large as few hundred lattice sites Q. In Fig. 1 


we present the phase diagram for the Hamiltonian (1) 
obtained by the QMC simulations. The phase boundary 
between the SF and Bose glass (BG) regions was found 
from the condition /i” 1 = tt(kPs) 1 ^ 2 = 3/2, wliere n is 
compressibility, p s is the superfluid stiffness. The criti- 
cal value for the exponent K at the SF-BG transition is 
known []:/[] to be 2/3. We see that the actual phase di- 
agram has little in common with that found in Ref. [Q, 
except at weak disorder where the phase boundary SF- 
BG was incorrectly identified with the SF-MI transition. 
To this end we note that simple considerations Q predict 
the upper bound for the MI-BG transition at A = E gap /A 
where E gap is the energy gap in a pure system. The cor- 
responding line is shown in Fig.l by open circles utilizing 
the data for E gap found in Ref. || . As expected, it is well 
below the SF-BG line [|j . We have explicitly checked that 
the State above the solid line at U = 2.2 is a gapless State 
with almost zero superfluid stiffness. 

In our opinion, the problem with the DMRG rnethod 
in disordered case arises because it implies that density 
matrix of any sufficiently large system block of R sites 
(. R » 1) is “representative”, i.e., describes atypical sys- 
tem behavior on the scale R; the variation of boundary 
conditions is supposed to be small when slightly increas- 
ing the block size. This assumption flatly contradicts 
the physics of strongly disordered case, the MI-BG tran- 
sition near the critical point, and the SF-BG transition 
in weakly interacting system. In the last two cases the 
large-scale behavior is governed by exponentially rare Lif- 
shitz regions of strong disorder, where the system wave- 
function changes drastically on rather small length scale, 
while at strong disorder such variations of the wave func- 
tion are just typical. DMRG rnethod, working with finite 
number of basis States fails to accomodate the sudden 
change in the boundary conditions when the block size is 
slightly increased, simply because the number of States 
kept is too small to reconstruct the correct new basis set. 
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